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We briefly review the growing efforts to set up a unified framework for the descrip¬ 
tion of neutrino interactions with atomic nuclei and nuclear matter, applicable in the 
broad kinematical region corresponding to neutrino energies ranging between few MeV 
and few GeV. The emerging picture suggests that the formalism of nuclear many-body 
theory can be exploited to obtain the neutrino-nucleus cross sections needed for both 
the interpretation of oscillation signals and simulation of neutrino transport in compact 
stars. 
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1. Introduction 

The description of neutrino interactions with nuclei and nuclear matter over a 
broad kinematical domain, besides being highly valuable in its own right, will be 
required as an input for the ongoing studies aimed at addressing two outstanding 
and fundamental physics issues: the violation of CP symmetry in the leptonic sector 
and the mechanism leading to supernovae explosions. 

Interactions of neutrinos of energy ranging from several hundreds MeV to few 
GeV determine the signals detected by many experimental searches of neutrino 
oscillations, 1 while the response of uniform nuclear matter to weak interactions 
at much lower energies, of the order of few MeV, plays a critical role in simula¬ 
tions of neutrino transport in compact stars. 2,3 In the region of high momentum 
transfer—typically |q| > 500 MeV—in which nuclear interactions can be described 
within the impulse approximation (IA) scheme, 4 the main difficulties involved in 


1 
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theoretical calculations of the neutrino-nucleus cross section arise from the neces¬ 
sity of combining a realistic model of nuclear structure and dynamics with a proper 
treatment of relativistic effects, which are known to be large. In the region of low 
|q|, on the other hand, the non relativistic approximation can be safely applied, 
but the scattering process may involve more than one nucleon, thus leading to the 
appearance of collective nuclear excitations. 

Electron scattering studies have provided ample evidence that in the IA regime 
the approach based on the factorisation ansatz and the spectral function formalism 
provides a viable computational framework to carry out accurate calculations of the 
nuclear inclusive cross sections. Of critical importance, in this context, is the use 
a dynamical model based on a realistic nuclear hamiltonian, allowing to take into 
account short range nucleon-nucleon (NN) correlations. Recently, the authors of 
Refs. 5, 6 have shown that, within the approach based on correlated basis function 
and the cluster expansion technique, the same hamiltonian can also be employed to 
obtain an effective interaction suitable for describing long range correlations—which 
are known to become important in the low energy regime—in a fully consistent 
fashion. 

The structure of the neutrino-nucleus cross section is outlined in Section 2, while 
Section 3 is devoted to a brief review of the dynamical model underlying Nuclear 
Many-Body Theory (NMBT). The formalism employed to describe the kinematical 
regimes corresponding to low and high energies are discussed in Sections 4 and 
5, respectively, while in Section 6 we summarise the main results and state the 
conclusions. 


2. The lepton-nucleus cross section 

For definiteness, we shall consider charged-current neutrino-nucleus interactions. 
However, the formalism outlined in this section can be readily generalised to the 
case of neutral current interactions. 6,7 

The double differential cross section of the process 


vg A —^ £ -f- X , 


(1) 


where A and X denote the target nucleus in its ground state and the undetected 
nuclear final state, respectively, can be written in the form (see, e.g., Ref. 8) 


d 2 a 

dk'dkQ 


G 2 F V 2 d |k'| 

16 7T 2 Ikl 


( 2 ) 


In the above equation, k = ( fc 0 , k) and k' = (k' Q , k') are the four momenta carried by 
the incoming neutrino and the outgoing charged lepton, respectively, k' = k'/|k' |, 
G F is the Fermi coupling constant and V ud is the CKM matrix element coupling u 
and d quarks. 

The tensor L^ is completely determined by lepton kinematics, whereas the 
nuclear tensor , containing all the information on strong interaction dynamics, 
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describes the response of the target nucleus to weak interactions. Its definition 

WV = £ (0|V|X) (X\J u \0) 5 {i \p 0 + q-p x ) , (3) 

x 

with q = (u, q) = k — k ', involves the initial and final states |0) and |X), with four 
momenta po and px , respectively, as well as the nuclear current operator, JA 
Note that the target ground state state |0) is independent of momentum transfer, 
while the state |X) includes at least one particle carrying momentum ~ q, and 
the current operator depends explicitly on q. As a consequence, a fully consistent 
theoretical calculation of the response tensor is only possible in the kinematical 
regime corresponding to |q|/m -C 1, with m being the nucleon mass, where the non 
relativistic approximation underlying NMBT is applicable. 9, 10 

On the other hand, the treatment of the region of high momentum transfer, 
relevant to event analysis of many neutrino experiments, requires a theoretical 
approach in which the accurate description of the nuclear ground state provided by 
NMBT is combined with a relativistically consistent description of both the final 
state and the nuclear current. 


3. Modelling nuclear structure and dynamics 

Nuclear Many-Body Theory (NMBT) is based on the tenet that nucleons can be 
treated as point like non relativistic particles, the dynamics of which are described 
by the hamiltonian 

fl = Elt+X>« + E < 4 > 

i—1 j>i =1 k>j>i =1 

In the above equation, p; is the momentum of the j-th nucleon, while the potentials 
Vij and Vijk describe two- and three-nucleon interactions, respectively. 

The nucleon-nucleon (NN) potential is obtained from an accurate fit to the 
available data on the two-nucleon system, in both bound and scattering states, and 
reduces to the Yukawa one-pion-exchange potential at large distances. State-of-the- 
art parametrizations of are written in the form 11 

18 

% = , ( 5 ) 

n= 1 

with rtj = |rj — Tj\ and 

Oif 6 = [1, (fTi ■ a-j), Sij\ <g> [1, (Ti ■ Tj)\ , (6) 

where er,; and r, are Pauli matrices acting in spin and isospin space, respectively, 
and 

3 

&ij ~ Z/r( (Ti ‘ r ij)( eT j ’ r ij) ~ ( <T, i <T i) ■ 
r ij 


(7) 
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The operators corresponding to n = 7,..., 14 are associated with the non static 
components of the NN interaction, while those corresponding to p = 15,..., 18 
account for small violations of charge symmetry. Being fit to the full Nijmegen 
phase-shift database, as well as to low energy scattering parameters and deuteron 
properties, the Argonne r>ig potential provides an accurate description of the two- 
nucleon system by construction. 

The inclusion of the additional three-body term, Vijk, is needed to explain the 
binding energies of the three-nucleon systems. 12 The derivation of Vijk was first dis¬ 
cussed in the pioneering work of Fujita and Miyazawa, 13 who argued that its main 
contribution originates from the two-pion exchange process in which a NN interac¬ 
tion leads to the excitation of one of the participating particles to a A resonance, 
which then decays in the aftermath of the interaction with a third nucleon. 

Commonly used models of the three-nucleon potential are written in the form 

Vijk = V*Z + Vg k , (8) 

where V) 2 £ is the attractive Fujita-Miyazawa term, while Vff k is a purely phe¬ 
nomenological repulsive term. The parameters entering the definition of the above 
potential are adjusted in such a way as to reproduce the ground state energy of 
the three-nucleon systems. Note that for A=3 the Schrodinger equation can still be 
solved exactly, using both deterministic or stochastic methods. 

The nuclear current consists of one- and two-nucleon contributions, the latter 
arising from processes in which the interaction with the beam particle involves a 
meson exchanged between the target nucleons. It can be conveniently written in 
the form 

^ = ^ + ^ = E^+E^- (9) 

* j>i 

The connection between the above current operator and the nuclear hamiltonian 
will be discussed in the next section. 

4. Nonrelativistic regime 

In the nonrelativistic regime, typically corresponding to |q| < 500 MeV, both the 
initial and the final state appearing in Eq. (3) are eigentstates of the nonrelativistic 
many-body hamiltonian of Eq. (4) , satisfying the Schrodinger equations 

ff|0)=£ o |0} , H\X)=E x \X). (10) 

In the case of light nuclear targets, the ground state wave function can be ob¬ 
tained from accurate stochastic approaches, such as Green’s Function Monte Carlo 
(GFMC). 14,15 In addition, the nuclear cross section can be conveniently rewritten in 
terms of the response functions ), obtained from Eq.(3) replacing the com¬ 

ponents of the current operator with their expressions obtained in the relativistic 
limit, taking into account terms up to order (|q|/m) 2 . 
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A fundamental feature of the description of neutrino-nucleus interactions at low 
and moderate momentum transfer is the possibility of employing a set of electroweak 
charge and current operators consistent with the hamiltonian of Eq.(4). 

The nuclear electromagnetic current, = (J° m , J em ), trivially related to the 
vector component of the weak current, is constrained by H through the continuity 
equation 16 


V ■ J em + i[H, J° m ] = 0 . 


( 11 ) 


Note that, because the NN potential does not commute with the charge operator 
J° m , the above equation implies that involves two-nucleon contributions, as 
shown in Eq.(9). 

The one-body electroweak operator is obtained from a non relativistic expansion 
of the covariant single-nucleon currents. Two-body charge and current operators 
are derived within the conventional meson-exchange formalism 1 ' ’ 18 or within the 
Effective Field Theory approach inspired to chiral perturbation theory (y;EFT). 19-23 
In this short review we will mainly discuss the former, in which the dominant 
static part of the realistic two-nucleon potential arises from exchange of effective 
pseudoscalar (7r-like) and vector (p-like) mesons, and the corresponding charge and 
current operators are projected out of the static components of the potential. As a 
consequence, the resulting vector current is conserved by construction. 

The conventional electroweak charge and current operators have no free parame¬ 
ters, except the nucleon-to-A axial coupling constant, which is fixed by reproducing 
the tritium Garnow-Teller transition strength in calculations based on the realistic 
hamiltonian discussed above. 

Nonrelativistic meson-exchange currents (MEC) have been used in analyses of 
a variety of electromagnetic moments and electroweak transitions of s- and p-shell 
nuclei at low and intermediate values of energy and momentum transfers. Taking 
MEC into account, a good agreement has been achieved between theoretical predic¬ 
tions and experimental data for the M 1 and E2 radiative transition rates between 
low-lying states, 24 ’ 25 /3-decays and electron- and muon-capture rates, 26,2 ' elastic 
and inelastic form factors measured in (e,e') scattering 10, 28-30 and radiative and 
weak capture reactions at low energies. 1 ' 


4.1. Quantum Monte Carlo approach 

Quantum Monte Carlo (QMC) methods were first applied to the study of proper¬ 
ties of light nuclei over three decades ago 31 (for a recent review of Quantum Monte 
Carlo methods for nuclear physics see, e.g., Ref. 32). Within the limits of applica¬ 
bility of NMBT, they provide a truly ab initio approach, allowing to perform exact 
calculations of a number of nuclear properties. 
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4.1.1. QMC calculations of the nuclear ground state 


The first calculations employed the variational Monte Carlo (VMC) technique, in 
which the stochastic Metropolis algorithm is used for evaluating the expectation 
value of a given many-body operator using a suitably parametrized trial wave func¬ 
tion, ’I'tv The parameters entering the definition of 'Ft are optimized by minimizing 
the variational energy 


Ey 


(^ T |iL|T r ) 
('F T |’F r ) _ 0 


( 12 ) 


which provides an upper bound to the ground-state energy Eq. It is worth noting 
that Monte Carlo methods can also be used in the search for the best variational 
parameters. 

Designing an accurate variational wave function requires a deep understanding 
of both the structure and dynamics of the nuclear system under consideration. 
Standard VMC calculations for light nuclei use a variational wave function of the 
form 


|^ T )=J-|$). (13) 

The long-range behaviour is described by the Slater determinant |<F). For example, 
in the case of uniform nucleon matter in the normal (i.e. non superfluid) phase, | $) 
is the wave function describing a non interacting Fermi gas. For light nuclei, on the 
other hand, |$) is usually written as a sum of Slater determinants, such as those 
resulting from from a small scale shell-model calculation. 

The short-range components of the wave functions are controlled by the corre¬ 
lation operator F , the structure of which reflects the complexity of the two- and 
three-nucleon potentials appearing in the nuclear hamiltonian 

i<j<k i<j 

In the above equation, S is the symmetrisation operator, needed to fulfill the re¬ 
quirement of antisymmetrisation of the trial wave function, while the two-body 
correlation operator exhibits a spin-isospin structure similar to that of the NN 
potential [compare to Eq.(5)] 

= ( 15 ) 

n 

implying that [Fij,Fik\ ^ 0. The scalar, f c = f 1 , and operator, / n>1 , pair cor¬ 
relation functions reflect the influence of the the short-distance behavior of the 
two-body potential and, at the same time, satisfy the boundary conditions implied 
by the requirement of cluster separability. 

As discussed below (see Section 4.2), reasonably accurate correlation functions 
are generated by minimising the two-body cluster contribution to the energy per 
particle. This procedure results in the derivation of eight Euler-Lagrange differential 
equations, involving a set of variational parameters. 33 
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Three-body correlation functions are induced by both the two- and three-body 
potentials. As for the latter case, the form suggested by perturbation theory 

Fijk = 1 + e x Vjjk (y rjj ,yr jk ,yrik), (16) 

X 

is usually employed. In the above equation, the subscript x labels the various con¬ 
tributions to the three-body force, the e x are small negative strength parameters, 
and y is a scaling factor. 

The GFMC method 14, 15 overcomes the limitations of the variational wave- 
function by using an imaginary-time projection technique to enhance the ground- 
state component of the starting trial wave function. The method relies on the ob¬ 
servation that 'P'r can be expanded in the complete set of eigenstates of the the 
hamiltonian according to 

l^r) = n) , H\n) = E n \n) , (17) 


which implies 


lim e- {H - Eo)T \* T ) =co|0>, (18) 

T—>-00 

where r is the imaginary time. Hence, GFMC projects out the exact lowest-energy 
state, provided Ty it is not orthogonal to the true ground state, i.e. Co ^ 0. 

The direct calculation of exp[— (H — Eq)t\ for strongly-interacting systems in¬ 
volves prohibitive difficulties . To circumvent this problem, the imaginary-time evo¬ 
lution is broken into N small imaginary-time steps, and complete sets of states are 
inserted, in such a way that only the calculation of the short-time propagator is 
required. This procedure yields the expression 

(R,v + 1 |e- ( "- £ ° )T |Ri> = J (Rv+i|e - ^ -Bo)Ar |Rjv)(RAr|e -( ' H_B ° )AT |Rjv-i) ■ ■ ■ 

x (R 2 \e- {H - Eo) \K 1 )dR 2 ...dR N , (19) 

where, for the sake of simplicity, the dependence on the spin-isospin degrees of 
freedom have been omitted. 

Monte Carlo techniques are used to sample the paths R, in the propagation. 
Note that, while being exact only in the At — > 0 limit, the accuracy of Eq.(19) can 
be checked by performing several simulations with smaller time step and extrapo¬ 
lating to zero. 

Because nuclear interactions are strongly spin-isospin dependent, the trial wave 
function is written as a a sum of complex amplitudes for each spin-isospin state of 
the system 


\^t) = JJi/>t(R) IXs(c)) | Xt(r)) ■ 


(20) 
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Within standard GFMC for nuclear physics applications, 34,35 the 2 A many-body 
spin states, defined as 

Xs=l I -i-l 3 -J^2 3 * * * 3 'I'A) 

Xs—2 = | 1 13 -i-2, • ■ • 3 !^) 

Xa=2* = | tl,t2,-..,t4) (21) 

are considered. The corresponding many-body isospin states can be obtained by 
replacing f and j, with p and n. Exploiting charge conservation, the 2 A isospin 
states can be reduced to A\/(N\Z\) states and, by assuming that the total isospin 
T is a good quantum number for the nucleus, the size of |xt( r )) can be further 
decreased. 

Because the GFMC imaginary-time evolution of Eq. (19) involves a sum over 
spin and isospin states at each step, the computing time grows exponentially with 
the number of particles. The largest calculations to date have been performed for 
the nucleus of 12 C and for the systems of 16 neutrons, corresponding to 540,672 
and 65,536 spin-isospin states, respectively. 

Over the past decade, the Auxiliary Field Diffusion Monte Carlo (AFDMC) 
has emerged as a more efficient algorithm for dealing with larger nuclear system. 36 
Within AFDMC, the spin-isospin degrees of freedom are described by single-particle 
spinors, the amplitudes of which are sampled using Monte Carlo techniques, and 
the coordinate-space diffusion in GFMC is extended to include diffusion in spin and 
isospin spaces. 

The early applications of AFDMC were based on variational wave function 
containing purely central correlation functions, and the accuracy of this approach 
for systems other than pure neutron matter was limited. Recently, the authors 
of Ref. 37 were able to add operator correlations to the trial wave function, and 
developed a novel importance sampling technique, making the accuracy of AFDMC 
comparable to that of GFMC even in systems containing both protons and neutrons. 

4.1.2. QMC studies of the response functions 

The calculation of the response functions involves major difficulties even in the 
region of |q| < 0.5 GeV and oj corresponding to quasi-free kinematics, where the 
consequences of the nucleon’s internal structure on nuclear dynamics can be sub¬ 
sumed into effective many-body potentials and currents. 


Integral properties of the response functions can be studied exploiting their sum 
rules, which are obtained from ground-state expectation values of appropriate com¬ 
binations of the current operators, thus avoiding the calculation of the full exci¬ 
tation spectrum of the target nucleus. GFMC calculations of the electromagnetic 
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Fig. 1. Euclidean electromagnetic response function of 12 C in the transverse channel at |q| = 570 
MeV. Experimental data are taken from Ref. 38. 


sum rules, 39 have demonstrated that a large fraction (~ 30%) of the strength in the 
transverse channel arises from processes involving two-body currents, and that in¬ 
terference effects between the matrix elements of one- and two-body currents play a 
major role. 40 These effects are typically only partially, or approximately, accounted 
for in existing perturbative or mean-field studies. 41-44 

The main drawback of the sum rules is that they do not provide any information 
on the distribution of strength; whether, for example, the calculated excess strength 
induced by two-body currents is mostly at large ui, well beyond the quasi-elastic 
peak, or it is also found in the quasi-elastic region. In addition, in the electro¬ 
magnetic case, comparison of theoretical and experimental sum rules is problem¬ 
atic, since longitudinal and transverse response functions obtained from Rosenbluth 
separation of the measured inclusive (e, e') cross sections are only available in the 
space-like region (w < |q|) and therefore must be extrapolated into the unobserved 
time-like region (w > |q|) before “experimental” values for the sum rules can be 
determined. 39 


Valuable information on the w dependence of the response functions can be inferred 
from the their Laplace transforms, also referred to as Euclidean responses, 40 defined 
as 

/»oo 

E^(q,T) = C fJ , u {q) dio e~ TU1 R^ v (q,uj) , (22) 

where Wth is the inelastic threshold and the C^ are ^-dependent normalization 
factors. In the case of the electromagnetic longitudinal ( L , or = 00) and 
transverse (T, or /j,v = xx) response functions, the normalization factors are 45 
Cl = Ct = 1/ \C p e {Q\ 0 )\ , where G P E is the proton electric form factor and 
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Fig. 2. Euclidean neutral-weak E xx (top panel) and E xy (lower panel) response functions of 12 C 
at |q| = 570 MeV. 


Qqe = IqI 2 — w qe) w qe being the energy transfer corresponding to quasi free kine¬ 
matics. In the neutral-weak response functions the normalisation factors are the 
same as those adopted in calculations of the sum rules. 46 

The chief advantage of the Euclidean response is that it can be expressed as a 
ground-state expectation value 

MNd) _ (0|Q£(q)e-( g - g °> T Off(q)|0) 

(0|e-(^-^o)-|0) ’ 1 j 

where H is the nuclear Hamiltonian, and Eq is a trial energy controlling the nor¬ 
malisation. 

In Fig. 1, taken from Ref. 47, the electromagnetic transverse Euclidean response 
function of 12 C, Et , is compared to the one obtained from the analysis of the world 
data carried out by Jourdan, 38 represented by the shaded band. The procedure fol¬ 
lowed to obtain the experimental Euclidean response is discussed in Ref. 47. Note 
that, in order to emphasise the large r behavior, the scaled Euclidean response 
-E M „(|q|,T) = exp[r q 2 /(2m)]i? /i ^(|qj,r) is displayed. The results obtained by in¬ 
cluding only one-body or both one- and two-body terms in the electromagnetic 
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transition operators are represented by open and solid circles, respectively. Two- 
body current contributions substantially increase the Euclidean response over the 
whole range of imaginary-time, thus implying that excess transverse strength is 
generated by two-body currents not only at ui > oj qe , but also in the quasi-elastic 
and threshold regions. The full predictions obtained including two-body currents 
are in excellent agreement with data. 

A similar enhancement brought about by the two-body currents has been ob¬ 
served in the Euclidean neutral-weak response functions, 47 displayed in Fig. 2. The 
neutral-current response i? X y(|q|,r) is due to the interference between the vector 
(VNC) and axial-vector (ANC) terms of the neutral current (NC), and in the inclu¬ 
sive cross section the corresponding -Ra^dqljW) enters with opposite sign depending 
on whether the process A(yi, v[) or A (A;, v[) is considered. 48 Hence, the difference 
between the neutrino and antineutrino cross sections turns out to be proportional 
to R xy . This difference may well have an impact on the determination of the CP- 
violating phase from neutrino and antineutrino events detected at DUNE. 49 

On the other hand, since for E xx (\q\,r) the interference between vector and 
axial-vector terms vanishes, the response is simply given by the sum of the terms 
with both transition operators arising from either the VNC or the ANC. For 
Exx{\q\, T ) these individual contributions, along with their sum, are displayed sep¬ 
arately. Both the E xx (\q\,r) and E X y(\q\,r) response functions obtained retaining 
one-body terms only in the NC are substantially increased when two-body terms 
are also included. This enhancement is found not only at low r, thus corroborating 
the sum-rule predictions of Ref. 46, but in fact extends over the whole r region. 
Moreover, the individual (VNC-VNC) and (ANC-ANC) contributions are about 
equally affected. 
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Fig. 3. (Color online) Transverse electromagnetic response functions of 4 He at |q| = 500 MeV. 
Experimental data are from Ref. 45. 
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The inversion of the Laplace transform, needed to retrieve the energy depen¬ 
dence of the responses, is long known to involve severe difficulties. A groundbreaking 
result has been recently reported by the authors of Ref. 47, who exploited the max¬ 
imum entropy technique to obtain the electromagnetic longitudinal and transverse 
responses of 4 He. 

Figure 3 shows the breakdown of the transverse response of 4 He at 
|q | = 500 MeV into one-nucleon current, two-nucleon current and interference 
contributions. 50 Note that the quantity displayed in the figure is normalized di¬ 
viding by the squared proton form factor. It clearly appears that including the 
two-nucleon currents leads to a sizable enhancement of the response, and that the 
large positive contribution of the interference term peaks at energy loss u> < uj qe . 
The agreement between the GFMC results and the data of turns out to be remark¬ 
ably good. 


4.2. Correlated Basis Functions and Custer Expansion Formalism 


The Green’s Function Monte Carlo method, while providing a most powerful com¬ 
putational scheme to carry out exact calculations of a variety of nuclear properties 
and scattering observables, is still limited to light nuclei, with A < 12. An alterna¬ 
tive approach, which has been extensively employed to study both medium-heavy 
nuclei 51 and nuclear matter, 52 is based on the use of correlated basis functions 
(CBF) and the cluster expansion technique (see, e.g., Refs. 53, 54, 55). 

Let us consider, for simplicity, uniform and isospin symmetric nuclear matter. 
In the is case, the correlated states are obtained from the corresponding states of 
the non interacting Fermi gas, |IV)fg, through the transformation 


I N) = 


F\N) 


FG 


(24) 


fg (N\ F\N)pq ' 

The operator F, embodying the correlation structure induced by the NN interac¬ 
tion, is written in the form (compare to Eq.(15)) 


F = Sl[F ij • 

ij 


(25) 


The two-body correlation functions Fjj , whose operator structure reflects the com¬ 
plexity of the NN potential, is written in the form 

6 

/b = E/ n M°b’ ( 26 ) 

n— 1 


including the contributions associated with the operators O”^ 6 of Eq.(6). 

The explicit calculation of matrix elements of a many-body operator, such as 
the nuclear hamiltonian H , between correlated states involves prohibitive difficul¬ 
ties, because it requires integrations over the coordinates—as well as summations 
over the discrete degrees of freedom—of many nucleons. This problem can be cir¬ 
cumvented expanding the matrix elements in series, the terms of which represent 
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the contributions of subsystems (clusters) containing an increasing number (2, 3, 
..., A) of particles. 

Within the cluster expansion approach, the expectation value of the hamiltonian 
in the correlated nuclear matter ground state, the minimum of which provides the 
variational estimate of the corresponding energy, can be written as 54 

E v = (0|H|0) = T 0 + ^(A E v )n , (27) 

n 

where the first term in the right hand side is the energy of the non interacting Fermi 
gas. The contributions to the cluster expansion (27) can be represented by diagrams 
and classified according to their topological structures. Selected classes of diagrams 
can then be summed to all orders, solving a set of integral equations referred to as 
Fermi Hyper-Netted Chain (FHNC) equations, 54 to obtain an accurate estimate of 
E v 

The correlation functions f n {rij) are determined from the minimisation of the 
expectation value of the Hamiltonian. The functional minimisation of the two-body 
cluster contribution to the energy per particle, (A£V) 2 , leads to a set of six Euler- 
Lagrange equations, to be solved with proper constraints that force f c and /( ra>1 ) 
to “heal” to one and zero, respectively. This is most efficiently achieved through 
the boundary conditions 56 

f n (r>d n )=6 n i , (28) 

dr 

where the healing distances d n are treated as variational parameters, to be de¬ 
termined from the minimisation of Ey. Additional variational parameters are the 
quenching factors a n which simulate modifications of the nucleon-nucleon poten¬ 
tial, arising from the screening induced by the presence of the nuclear medium, and 
the set of scaling factors /3 n , often applied to the correlation functions /”. 

4.3. The CBF effective interaction 

The applications of the formalism of correlated basis function to the study of the 
nuclear matter response at low momentum transfer exploit the CBF effective inter¬ 
action 

= ( 29 ) 

j>i 

defined by the relation 57 

(H) = (0|ff|0> =T 0 + FG (0|E eff |0) FG . (30) 

The left hand side of the above equation, computed within the variational approach 
using the FHNC summation scheme, is assumed to provide a good approximation 
to the ground state energy, while the right hand side is computed at low order of the 
cluster expansion, i.e. two- or three-body cluster level. The range of the correlations 
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is adjusted in such a way as to satisfy Eq. (30), which implies that the low order 
calculation reproduces the variational ground state energy. 



p (fm- 3 ) 


Fig. 4. Density dependence of the energy per particle of isospin symmetric nuclear matter in the 
low density regime, computed with a bare hamiltonian containing the Argonne v' 6 potential and 
the UIX three-body interaction model. The solid line displays the results obtained using the CBF 
effective interaction, while the full variational energy per particle is represented by the shaded 
region, accounting for the uncertainty arising from the treatment of the kinetic energy term. 5 


The authors of Ref. 57, who performed the first calculation of V eS of Eq.(30), took 
into account two-body cluster diagrams only. This prescription, while leading to a 
very simple and transparent expression of V efc , fails to account for three-body forces, 
which are long known to play a critical role in determining the energy spectrum 
of light nuclei, as well as the saturation properties of isospin symmetric nuclear 
matter. 

The approach proposed in Ref. 58, in which interactions involving more than 
two nucleons are included through a density dependent modification of the NN 
potential at intermediate range, has been implemented in the effective interaction 
formalism in Refs. 59, 60. More recently, the CBF effective interaction has been 
substantially improved by the authors of Refs. 5, 6, who explicitly included the 
three-body cluster contributions, allowing for a more realistic description of three- 
nucleon forces at microscopic level. The energy-dependence of the binding energy 
per nucleon in isospin symmetric matter obtained from the effective interaction 
of Refs. 5, 6, which includes the effects of the UIX model of the three-nucleon 
potential, is illustrated in Fig. 4 (taken from Ref. 5). It turns out to be fairly close 
to that resulting from the full FHNC calculation, and exhibits saturation, albeit at 
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density larger than the empirical value po = 0.16 fm~ 3 . 


4.4. Weak response of nuclear matter at low momentum transfer 

The cooling of neutron stars is driven by the energy loss caused by the flux of 
neutrinos leaving the star. This effect can be conveniently parametrized in terms of 
the neutrino mean free path (NMFP), which is one of the critical inputs required 
for large-scale simulations of neutrino transport. 61 

Neutrino and antineutrino interactions in neutron matter are also relevant for 
understanding the evolution of the very neutron-rich matter formed in neutron- 
star mergers, since they can potentially affect the neutron to proton ratio and 
significantly impact the r-process in neutron star mergers, currently considered to 
be an important source for r-process nucleosynthesis. 

The relatively low momentum scale of the above processes, typically |q| < 50 
MeV, allows for a nonrelativistic treatment of both the current operators and the 
final states entering the hadronic tensor. Earlier calculations have accounted for 
the effects of correlations in the nuclear wave function through empirical effective 
interactions or Landau parameters, derived from Skyrme-like effective interactions, 
and using the Random Phase Approximation (RPA). 62,63 However, the effect of 
correlations on the current operators was totally neglected, as bare weak operator 
were considered. It is well established that such approach is inconsistent, and that 
a consistent set of effective operators and effective interactions must be included in 
a more accurate treatment of the nuclear response functions. 

The CBF effective interaction approach is best suited to define effective weak- 
current operators that are consistent with the effective interaction of Eq. (30). 
Under the assumption that the nonrelativistic final states entering the hadronic 
tensor can be described by CBF states of Eq. (24), the effective operators can be 
defined through their transition matrix elements 


Wf|o) 


FG (X|J-tJ M J-|0) FG 

\J FG(0|^ rt ^ r |0) FG FG (A|J-tJ-|A) FG 


FG {X | 10) FG . 


(31) 


Existing calculations based on the CBF effective interaction approach have only 
accounted for transitions between the correlated ground-state and correlated one 
particle-one hole (lplh) excited states, which amounts to setting FG (X| |0) FG ~ 
fg (ph\ J“ ff |0) FG , where p and h denote both the momentum and the spin and isospin 
projections specifying single nucleon state. 

The effective operators encompass short-range correlations, but fail to account 
for long range correlations, responsible for collective modes. The |p/i) FG states are 
not eigenstates of the effective hamiltonian and, as a consequence, there is a residual 
interaction that can induce transitions between different lplh states. 

The effect of long-range correlations can be included in the effective interaction 
formalism using the Tamm-Dancoff (TD) approximation, i.e. expanding the final 
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state in the basis of one lplh states according to 



Wtd=£ Cp h \ph) FG . 

(32) 


ph 


The excitation energy 

Ex of the state |X), as well as the coefficients 

Cp h , are 

determined by solving 

the eigenvalue equation 


H 0 s V)td = Ex |X) T d , 

(33) 

with 




"•*= -x:S+£<- 

(34) 


i j>i 


Note that the correlations defining the effective interaction used in the previous 
equations are the same appearing in the definitions of the effective current operators. 

The effect of long-range correlation is apparent in the spin-response of pure 
neutron matter, 6 displayed in Fig. 5, (taken from Ref. 6). The dashed and dot-dash 
lines represent the spin-transverse (pu = xx + yy) and spin-longitudinal (pv = zz) 
response functions of pure neutron matter at density p = 0.16 fm -3 computed in 
correlated Tamm-Dancoff (CTD) approximation, respectively. For comparison, the 
solid line corresponds to the spin-density response obtained from the Landau theory, 
the parameters of which have been consistently derived from the same effective 
interaction. 64 
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Fig. 5. Spin-transverse (dashed line) and spin-longitudinal (dot-dash line) responses of pure neu¬ 
tron matter, computed within the CTD and CHF approximations at p = 0.16 fm —3 for momentum 
transfer |q| = 0.1 fm —1 . The solid line has been obtained from Landau theory, according to the 
approach of Ref. 64 . 


When long-range correlations are accounted for, the peak associated with the 
collective excitation sticks out in both the spin-longitudinal and spin- transverse 
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channels. On the other hand, in the correlated Hatree-Fock (CHF) scheme, in which 
nuclear correlations only enter via the effective operators and the quasiparticle 
energies, the strength of the response function is smoothly distributed over the 
particle-hole continuum. 

Long-range correlations have been shown to produce a similar effect in the Fermi 
and Gamow-Teller responses of isospin-symmetric nuclear matter , 5,5 ‘ ’ 60 whereas no 
collective mode is observed in the neutron matter density response function. 1 ’ 

The neutrino mean free path for low energy neutrino scattering and neutrino 
absorption processes in cold isospin symmetric matter has been found to be largely 
affected by both short- and long-range correlations. 57 For densities ranging from 
Po/2 to 3/2po, the NMFP obtained from the CTD response functions is ~ 2.5 — 3.5 
times larger than the one of the noninteracting Fermi Gas case. In addition, the 
NMFP for scattering is 2 times larger than that for absorption, indicating that the 
cross section for charged-current transitions is 2 times larger than the one associated 
with neutral-current process. 

The role of long-range correlations in determining the NMFP associated with 
neutrino scattering processes in cold neutron matter has been investigated within 
Landau theory in Ref. 64. Figure 6, taken from Ref. 65 shows the density dependence 
of the mean free path of a non degenerate neutrino with an energy E = 1 MeV. 
The results of Landau theory, obtained including tensor interaction terms (solid 
line) and neglecting them (open circles), are compared with those corresponding to 
a free neutron gas (dot-dashed line). It clearly appears that inclusion of interaction 
effects leads to a large enhancement of the NMFP over the whole density range. 
The collective mode of the spin-density response function increases the scattering 
cross section, hence reducing the NMFP, by about 25%. It is worth nothing that 
the the results of Fig. 6 have been confirmed by the CTD calculations reported in 
Ref. 6. 

5. Relativistic regime 

The dynamical model discussed in Section 3 and the formalism of correlated basis 
functions can be also employed to describe the nuclear response in the kinematical 
region of large momentum transfer, in which |q| _1 -C d , with d being the average 
NN separation distance in the target. In this regime nuclear scattering can be rea¬ 
sonably assumed to reduce to the incoherent sum of elementary scattering processes 
involving individual nucleons, and the IA is expected to be applicable. 

Within the IA scenario, the difficulties associated with the relativistic treatment 
of the nuclear final state and current operator are circumvented exploiting the 
factorisation ansatz , which amounts to: (i) neglecting the contribution of the two- 
nucleon current, and (ii) rewriting the nuclear final state in the form 

\X) = |p) 0 |n A -i) Pn) • (35) 

In the above equation, the state |p) describes a non interacting nucleon, while 
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p [fm 3 ] 


Fig. 6. Density dependence of the mean free path of a nondegenerate neutrino with energy E = 1 
MeV in cold neutron matter. 64 Landau theory results obtained including (not including) tensor 
terms are represented by the solid line (open circles). For comparison, Fermi gas results are also 
shown, by the dot-dashed line. The dashed line refers to the NMFP obtained without accounting 
for the collective mode in the spin-density channel. 


|riA-i,Pn) is an eigenstates of the nuclear hamiltonian of Eq.(4), describing the 
recoiling (A — l)-nucleon system with momentum p„. 

The use of Eq. (35) allows to rewrite the nuclear transition matrix element 
in a most simple and transparent form, consisting of the matrix element of the 
one-nucleon currents between free nucleon states—which can be computed exactly, 
retaining the fully relativistic expressions of the currents—and the nuclear am¬ 
plitude involving the target ground state and the state of the recoiling spectator 
system—which can be safely obtained from non relativistic many-body theory, since 
no nucleons carrying large momenta, ~ q, are involved. 00 

The resulting expression of the differential nuclear cross section is 66 

dajA = J d 3 k dE Pi(k,E) dai , (36) 

i 

where dai is the corresponding cross section describing scattering on the i-th nu¬ 
cleon, the momentum and removal energy of which are distributed according to the 
spectral function P(ik,E). 

Accurate theoretical calculations of the nuclear spectral function have been car¬ 
ried out for the few nucleon systems, with A < 4, as well as for isospin symmetric 
nuclear matter (see Ref. 66 and references therein). For intermediate mass nuclei, 
realistic spectral functions have been constructed within the Local Density Approx¬ 
imation (LDA), in which the results of nuclear matter calculations are combined 
with the empirical information extracted from the measured (e, e'p) cross sections. 67 

It is very important to realise that, owing to the presence of strong nucleon- 
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nucleon correlations in the nuclear ground state, the recoiling (A—l)-nucleon system 
is not necessarily left in a bound, one hole, state. Two hole-one particle states, in 
which one of the spectator nucleons is excited to the continuum, typically contribute 
15-20% of the spectral function normalisation, the corresponding strength being 
located at large momentum (|k| > 300 MeV) and energy (E > 40 MeV), well 
outside the region corresponding to shell model states. 



i/ [GeV] 

Fig. 7. Top: inclusive electron scattering cross section oft carbon at beam energy E e = 1.3 GeV and 
electron scattering angle 9 e = 37.5 deg, corresponding to Q 1 = 0.4 GeV 2 at the A production peak, 
as a function of the electron energy loss v. The solid line corresponds to theoretical calculations 
carried out using the proton and neutron structure functions of Ref. 68, while the shaded region 
has been obtained with those resulting from the analysis of Ref. 69. The data are taken from Ref. 
70. Lower panel: same as in the upper panel, but for oxygen target, E e = 1.2 GeV and 9 e = 32 deg, 
corresponding to Q 2 = 0.26 GeV 2 at the A production peak. The data are taken from Refs. 71, 
72. 


5.1. Spectral function formalism 

To the extent to which the target spectral function is available, Eq. (36) can be used 
to perform theoretical calculations of the nuclear cross section within the IA. The 
other required input, i.e. the nucleon cross section dcq, can be obtained—at least 
in principle—from hydrogen and deuteron data. This procedure has been widely 
employed to analyse the large body of inclusive electron scattering data, both in 
the quasi elastic sector and beyond pion production threshold. As an example, 
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Figs. 7 and 8 show the inclusive cross section corresponding to different targets and 
kinematical setups. 

Inspection of Figs. 7 and 8 indicates that the approach based on the IA and the 
spectral function formalism provides a good description of the data in the quasi 
elastic sector, corresponding tow« wqe, in which the elementary electron-nucleon 
cross section can be written in terms of the proton and neutron vector form factors. 
On the other hand, the results of Fig. 7 suggest that the available parametrisations 
of the nucleon structure functions in the region in which the excitation of the 
A resonance is the dominant reaction mechanism involve a significant degree of 
uncertainty. At larger momentum transfer deep inelastic scattering clearly appears 
to take over at wqe, and the agreement between theory and data over the 
whole range of energy loss turns out to be remarkably good. 

5.2. Corrections to the impulse approximation 



CJ [GeV] 


Fig. 8. Inclusive electron scattering cross section off 3 He at beam energy E e = 11 GeV and electron 
scattering angle 9 e = 8 deg, corresponding to Q~ = 2.1 GeV 2 at the quasi elastic peak, as a func¬ 
tion of the electron energy loss u. The solid line shows the full theoretical result, while the dashed 
and dot-dashed lines correspond to the quasi elastic and inelastic contributions, respesctively. The 
data are taken from Ref. 73 (Adapted from Ref. 74). 


The capability of obtaining accurate estimates of the nuclear cross section in 
the quasi elastic channel over a broad kinematical range is important for the in¬ 
terpretation of the signals detected by many neutrino experiments. For example, 
quasi elastic scattering provides the dominant contribution to the event sample 
collected by the MiniBooNE Collaboration using a neutrino flux of mean energy 
~ 800 MeV. 75 
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240 MeV, 36 deg 200 MeV, 60 deg 240.4 MeV, 60 deg 







Fig. 9. (color online). Double differential electron-carbon cross sections in the QE channel of 
Ref. 76, compared to the data of Ref. 78, 79, 80. The solid lines correspond to the result of the 
full calculation, whereas the long-dashed lines have been obtained neglecting FSI. The difference 
between the solid and short-dashed lines illustrates the effect of using alternative treatments of 
Pauli blocking. For comparison, the predictions of the Fermi gas model are also shown, by the 
dotted lines. The panels are labeled according to beam energy, scattering angle, and values of |q| 
and Q 2 at the quasielastic peak. 


The comparison between the he results of Ref. 76 and the measured electron 
scattering cross sections (see Figure 9), suggests that a remarkably good agreement 
between theory and data in the region of the quasi elastic peak can be achieved 
correcting the IA results to take into account final state interactions (FSI) between 
the struck nucleon and the spectators, the effects of which can be described within 
a generalisation of the spectral function approach, discussed in Ref. 77. 
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Additional corrections to the IA scheme arise from processes involving two- 
nucleon currents. It has been suggested 42,81 that the inclusion of these processes, 
which have been shown to play an important role in the non relativistic regime (see 
the discussion of Section 4.1.2) may in fact explain the large disagreement between 
the results of Monte Carlo simulations and the double differential neutrino-carbon 
cross section measured by the MiniBooNE Collaboration.' 4 The generalisation of 
the factorisation ansatz proposed in Ref. 82, 40 will provide a consistent framework 
to carry out calculations of nuclear amplitudes involving two-nucleon currents, com¬ 
bining the fully relativistic expression of the currents and a description of nuclear 
dynamics taking into account short range correlations. 


6. Summary and perspectives 

The results of extensive studies carried out over the past decade provide convincing 
evidence that the model of nuclear structure and dynamics based on non relativistic 
many-body theory can be exploited to achieve a consistent description of neutrino 
interactions with nuclei -and, more generally, nuclear matter—over an energy range 
spanning three orders of magnitude. The availability of such a description will 
be needed to reach the level of accuracy required by the planned experimental 
searches of CP violation in the leptonic sector, as well as to significantly improve 
the modelling of neutrino transport involved in large scale simulations of compact 
star evolution. 

The Monte Carlo approach, while being capable to yield nearly exact results, 
is—at least for the foreseeable future—limited to the quasi elastic sector and the 
non relativistic regime. Moreover, its extension to nuclei heavier than carbon would 
require the use of enormous computing resources. On the other hand, approxima¬ 
tion schemes based on the same dynamical model, described by a realistic nuclear 
hamiltonian, provide a viable alternative approach, suitable for treating the regime 
of high momentum and energy transfer, in which relativistic effects are important 
and the hadronic final states involve hadrons other that protons and neutrons. 

The picture emerging from the available results suggests that Monte Carlo 
techniques may be effectively combined with more approximated approaches, to 
both improve their accuracy and widen their scope. As an example, Monte Carlo 
estimates—even at variational level—of the nuclear amplitudes involved in spectral 
function calculations may provide valuable information, useful to reduce the theo¬ 
retical uncertainty associated with the use of the local density approximation. In 
addition, the equation of state of cold nucleon matter computed using the Monte 
Carlo approach will provide more precise benchmarks for the determination of the 
effective interactions within the formalism of correlated basis function. 

The development of a unified description of neutrino interactions based on a 
realistic model of nuclear interactions—strongly constrained by phenomenology— 
appears to be possible, and well on its way. 
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